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Abstract. The currently used linear and nonlinear lens 
inversion techniques are based on distortion estimators 
whose complicated source redshift dependence makes the 
influence of the redshift distribution of the sources diffi- 
cult to take into account and to analyze. However, the 
lens equations can be explicitly averaged over the source 
redshift distribution by a suitable choice of the estima- 
tors. Lens reconstruction procedures are outlined, which 
make use of all the information and in which all the un- 
known quantities of the problem can be recovered simulta- 
neously in both the linear and nonlinear regime, for either 
noncritical or critical lenses. These procedures require no 
prior knowledge of the redshift distribution of the sources, 
but are possibly dependent on the Cosmology in some in- 
stances. Possible methods of recovery of the source redshift 
distribution itself are briefly discussed at the same time. 
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k* - ■ 1. Introduction 



A decade after the discovery of giant gravitational arcs in 
distant clusters of galaxies (Soucail et al. 1987, Lynds and 
Petrosian 1986), the strong and weak lensing distortions 
induced by the distribution of foreground dark matter on 
the position and shape of background sources have become 
one of the major tools in modern observational Cosmology. 
Early work (e.g., Kochanek 1990, or Miralda-Escude 1991, 
1993) has focused on parametric reconstruction of the lens 
properties; a significant progress along this line has been 
performed by the identification and use of multiple im- 
ages of a same source (e.g., Mellier et al. 1993, Kneib et 
al. 1993, 1995). Two major breakthroughs have been made 
in the recent years: on the theoretical side, powerful linear 
and nonlinear nonparametric lens reconstruction methods 
have been proposed (e.g., Kaiser and Squires, 1993; Kaiser 
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1995; Seitz and Schneider 1995), while on the observa- 
tional side, new and impressive techniques have been de- 
veloped to measure the distortion field in the weak lensing 
regime from the distortion of individual sources (Bonnet 
et al. 1994) and from the autocorrelation function of whole 
images (Van Waerbeke et al. 1997). Finally the end of the 
long cosmological quest for the magic numbers of the Uni- 
verse (f2 and A; see, e.g., Fort et al. 1997) and even the 
determination of the power spectrum of large scale struc- 
tures (Bernardeau et al. 1997) seem within reach of the 
lensing approach. 

However, a remaining vexing difficulty plagues these 
attempts, namely the complicated dependence of the dis- 
tortion estimators on the background source redshifts, 
which makes the influence of the source redshift distri- 
bution on the lensing analyses difficult to quantify (see, 
e.g., Luppino and Kaiser, 1996). The purpose of this Let- 
ter is to show how this difficulty is resolved by an ap- 
propriate choice of the distortion estimators (this was al- 
ready partially pointed out by Kochanek 1990; however, 
the implications for the elaboration of reconstruction pro- 
cedures with a distribution of source redshifts seem to have 
gone unnoticed; for an alternative approach, see Seitz and 
Schneider 1997). Sect. 2 is devoted to an exposition of the 
lensing and deformation equations used in the paper, both 
to specify the adopted notation and introduce the distor- 
tion estimators. Sect. 3 begins with a presentation of the 
relevant probability distributions of images and sources, 
and moves on to the statistical properties of the chosen 
distortion estimators. Sect. 4 discusses these results and 
concludes this Letter by outlining possible methods of re- 
covery of the lens projected mass distribution and of the 
source redshift distribution for noncritical as well as crit- 
ical lenses. 



2. Lensing equations and notations: 

In what follows, we make use of quantities which can be 
defined in similar ways for the sources, the images and 
the lens. Such quantities are accordingly indexed with a 
subscript or a superscript "I", "5"' and "L". Latin sub- 
scripts refer to the angular directions in the plane of the 
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sky, which are labeled x\ and xi- With these definitions, 
the lens equation reads 



x s = x 7 - V^x'), 



(1) 



where the bidimensional lensing potential 4> is related to 
the lens surface density £ by the bidimensional Poisson 
equation 



(2) 



The critical surface density S c just introduced is defined 
by 



2irG D ls D l ' 



(3) 



where D L , D s , and D LS are the angular-diameter dis- 
tance of the observer to the lens, the observer to the 
source, and the lens to the source, respectively. 

Extended objects are characterized by their surface 
brightness Is{x-s) = -M x/ )> tne position of their centers 
in the source plane and in the image plane 



JJ Ii,s{^' S ) d 2 * 1 ' 3 ' 
and their shape matrix 

Ml , s II iiA*i s -x!Vs s -*!•*) d 2 ^ 3 



(4) 



(5) 



Following Kneib et al. (1994), it is convenient to ex- 
press the normalised source and image shape matrices 
-My(detM) 1 / 2 in terms of the orientation of the principal 
axes, 9 (defined such that tan(20) = 2M 12 /(M n - M 22 ) ), 
and the major axis a and minor axis b of the equivalent 
ellipse as 



M 



5 + |t| cos(26») |r|sin(26>) 



(detM) 1 ^ ^ |r|sin(20) 5 - |r| cos(20) J ' 



(6) 



where the modulus of the complex quantity r = |r| exp 2i9 
is given by 



2a b 



(7) 



and 

5= ■ 



t T "0 



2a 6o 



(1 + M 2 ) 1/2 . 



Finally, one can define a last observable, a, which mea- 
sures the amplitude of the deformation matrix: 



a = 8(detM)^ 2 = 2ab. 



(9) 



When the distortions are moderate (so that the amplifica- 
tion matrix can be assumed to be nearly constant across 



the image), the source and image shape matrices M are 
related by 

M s = a^M'a- 1 , (10) 
where the inverse of the amplification matrix a -1 reads 



1-0, 

-0,12 



11 -<P,12 
1 - 0,22 



(11) 



The coma designates partial derivatives with respect to 
coordinates represented by their indices. 

It is worth pointing out here that Van Waerbeke et 
al. (1997) have introduced a different distortion estima- 
tor, the second moments of the autocorrelation function 
in the image plane. This new matrix estimator possesses 
at least two fundamental advantages over the more con- 
ventional individual image shape matrix: first, it makes 
maximal use of all the information contained in CCD im- 
ages, including the information contained in the noise, so 
that this autocorrelation matrix can be extracted from 
the data everywhere in the field with unprecedented ac- 
curacy, which makes it the most promising data analysis 
tool to this date; furthermore, this estimator transforms 
like the shape distortion matrix introduced above, so that 
the analysis developed below can be literally transposed 
to it. 

Introducing the convergence 



k = 



A(f> 



S 



and the complex shear 
1 



7 = |7|exp(i2# L ) = ^(0, 2 2 



,ll) - 10,12, 

the inverse amplification matrix reads 



(12) 



(13) 



„-i_ A-«+|7|cos(20i) |7|sin(20 L ) \ 
|7|sin(20 L ) l-n-\ 7 \cos(29 L ) J ■ 

Besides 7, two other quantities related to the lens poten- 
tial, tl and 5l, are needed in the following analysis, and 
write 



TL = 



-2(1-«)| 7 | 



-2(1 -«)N 



\A~i\ |(1_«)2_| 7 |2|' 

(i- K ) 2 + bf (i-«o 2 + | 7 | 2 



L4-i| 



:i-«) 2 -mt 



(15) 
(16) 



^ where A 1 = det(a 1 ). Note that 8\ = 1 + r|, as the 



equivalent source and image quantities. 

With these definitions, the shape matrix transforma- 
tion equation [Eq. (10)] yields 



07 = 



(17) 



\T S \sm(29 s -29 L ) = \ Tl \ sin(20 7 - 29 L ), (18) 
|r s I cos(2# s - 29 L ) = S L \ TI \ cos(2#/ - 29 L ) - t l Si. (19) 
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The last equation can be inverted, which gives 

M cos(20/ - 26 L ) = 5 l \t s \ cos(2# s - 28 L ) + t l S s , (20) 

In the rest of the paper, the complex distortion parameter 
t and the distortion amplitude a are used instead of the 
shape matrix. 

These distortion equations present several advantages 
over the ones which are presently used in inversion pro- 
cedures. First, the linear relation between source and im- 
age distortion estimators simplifies greatly the averaging 
over the source orientation. Second, the \A~ X \ denomina- 
tor makes an explicit redshift averaging of the equations 
possible, as already pointed out by Kochanek (1990), and 
shown in Sect. 3. On the other hand, r is sensitive to the 
intrinsic ellipticity distribution of the sources, so that a 
single highly elongated source can in principle bias the es- 
timation of the local averages used in the next sections; 
however, this seems likely to happen only rarely in prac- 
tice for reasonable source redshift distributions. 



3. Probability distributions and distortion statis- 
tics: 

The lens and the distortion equations involve both deter- 
ministic paramaters (the lens redshift and mass distribu- 
tion, the cosmological parameters) and random variables 
(the position x s , the complex distortion ts, the distor- 
tion amplitude as and the redshift zs for the sources, the 
position x J , the complex distortion 77 and the distortion 
amplitude 07 for the images; the number of sources and 
images per unit area is also a random variable, but it is 
not used here). It is commonly assumed that the positions, 
distortion estimators, and redshifts of the sources (i.e., in 
the absence of the lens) are all independent random vari- 
ables (however, the following argument is unaffected by 
a possible correlation between 175 and \ts\). Furthermore, 
the positions and orientations 9s of the sources are as- 
sumed to be uniformly distributed random variables. 

With the assumption just recalled, the joint 
probability density distribution of source parameters 
p(-x s ,T S ,as,zs) = {^/^)p{ts)p{<js)p{zs)- The statisti- 
cal quantities of interest are the averages of the distortion 
estimators over the source random variables, at a given po- 
sition in the lens plane, i.e., the probability distribution 
needed is the conditional probability p(rs, <js, zs/xl 1 ) = 
pi* 1 ,T S ,<7s,zs)/vW) = p(x ! 7 zs)p{t s )p{(Ts)/p(x 1 )- Fur- 
thermore, 



p(x',^) = i-p(^)|^- 1 | 
from which one obtains 
p^ I ) = -L(\A- 1 \) zs , 



(21) 



where (|A _1 |) ZS stands for the redshift average of \A~ 1 \ at 
given image position x 7 . Finally, the required probability 
distribution reads 



P(rs, <?s, z s /x r ) = p(T S )p(as)p(z s ) 



(\A^\), 



(23) 



For critical lenses, J p(x 7 ) cPx 1 > 1, because multiple 
images are absolutely correlated events and should not be 
coadded as independent ones. If multiplied by the number 
of sources in the source plane, Eq. (22) yields the expec- 
tation value of the number density in the image plane. 

Averaging the distortion equations for Si, the com- 
plex distortion 77 and the distortion amplitude 07 over 
the sources at a given position in the image plane now 
yields 



(Si) 
(n) 



( 5 s)\t s \ 
( S s)\t s \ 



(Sl\A-i\) zs 

i\A-i\) zs ' 
(tl\A-'\) zs 7 



(24) 
(25) 
(26) 



where 77 = |r/| exp(2i#/), and where the ratio 7/I7I does 
not depend on the redshifts of the sources. The left-hand 
sides of Eqs. (24) through (26) are provided by the obser- 
vations. 

It is necessary to express the redshift averages appear- 
ing in the preceding equations in more detail. To this pur- 
pose, let us introduce the lens critical radius b = D^s/Ds, 
and define a redshift independent convergence K and 
shear T by 



K = bK, 

so that 

(T L \A-i\) zs = -2(b)\r\ + 2{b2)K\r\, 

(5 L \A- 1 \) zs =l-2(b)K + (b 2 )(K 2 + \T\ 2 ). 



(27) 



(28) 



(29) 



(30) 



The sign of A^ 1 depends on the position in the image 
plane and on the source redshift. If the sources are all at 
the same redshift, the sign change occurs on the critical 
lines, i.e. a zero- width region. If the sources have some dis- 
persion in redshift, the critical lines become critical zones. 
Two cases naturally emerge, depending on the size of these 
critical zones. 

If the critical zones are too small to be of practical 
importance, one has 



( 22 ) (lA-'Dzs = sgn^- 1 ) [1 - 2{b)K + (b 2 )(K 2 \T\ 2 )] ,(31) 
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where sign(A _1 ) is well-defined outside the (neglected) 
critical zones. Both (b) and (b 2 ) depend only on the dis- 
tribution of source reshifts and on the Cosmology. Also, 
{5 I ) 2 /{S s ) 2 = l + \(r I )\ 2 /{S s ) 2 . 

Otherwise, at every position x 7 one can define a red- 
shift domain Z + for which A^ 1 > 0, and a complemen- 
tary redshift domain Z~ where A^ 1 < 0; these domains 
are known if the lens surface density is known, for a given 
Cosmology. Let us define = J z± b n p(z$) dzs; one has 
B+ + B~ = {b n ) (= 1 for n = 0), where (b) and (b 2 ) 
have the same (position independent) value as in regions 
of constant sign of A~ x (e.g., outside the asymptotic ex- 
ternal critical line). Then 

(lA- 1 !)^ = (AB ) - 2(AB 1 )K + (AB 2 )(K 2 - |r| 2 ), (32) 

where AB n = £?+ — B~ . Note that the three quantities 
AB n are all related to each other through the source red- 
shift distribution and the Cosmology, so that the spatial 
dependence of ({A^ 1 \} zs provides us with an integral equa- 
tion for p(zs), for a given Cosmology and for a given lens 
surface density distribution. As cover the whole red- 
shift range of the distribution when one crosses the critical 
zones, this integral can in principle be solved for p(zs), 
possibly in a parametric way. 

4. Reconstruction procedure and discussion: 

Following a lead by Kaiser and Squires (1993), Seitz and 
Schneider (1995) have shown that the convergence and the 
complex shear are related through an invertible convolu- 
tion equation which yields 

r(x 7 ) = I J d 2 y J G(x 7 - y'My 7 ), (33) 
K^ 1 ) = I J d 2 y J Re (G(x 7 - y 7 )r*(y 7 )) , (34) 
where the kernel G is given by 

G(x) = x "-^ + 2 ^ 2 (35 ) 
l x l 

Two different cases must be distinguished, depending on 
whether one can ignore the critical zones defined above. 
Let us first assume that we can, which is probably cor- 
rect in a majority of cases. Then Equations (24) through 
(35) [especially Eqs. (24), (26) and (33)] provide us with 
enough independent constraints for the five quantities (K, 
Re(r), Im(r), (b) and (& 2 )), which can therefore in prin- 
ciple be solved for as functions of x 7 (with appropriate 
noise filtering; Kaiser and Squires, 1993, and Seitz and 
Schneider, 1995). Note that (|yl _1 |}^ s is obtained from the 
distortion amplitude a but any other estimator of \A~ X \ 
would do as well (e.g., it is possible that the magnification 
might be extracted with more accuracy from the image au- 
tocorrelation itself than from its second order moments). 
Also, a and S are expected to be observationally less well 



constrained than r. However, the expected constancy of 
(b) and (b 2 ) in the image plane provides us with an ad- 
ditional constraint on the inversion procedure, as will the 
identification of multiple images, if any, so that, although 
this inversion procedure would probably be difficult to 
implement from the observed distortions of individually 
identified objects alone, it seems within reach of the au- 
tocorrelation method proposed by Van Waerbeke et al. 
(1997). Note however that the strong amplification near 
the critical zones biases the observed distribution of ob- 
jects towards fainter or farther sources, which can possibly 
invalidate the assumed constancy of (b) and (b 2 ). 

Note also that the local degeneracy pointed out by 
Schneider and Seitz (1995) is broken as long as (b 2 ) is 
significantly different from (b) 2 . As expected, the global 
degeneracy of the distortion equations cannot be raised 
without the use of external information. In the same vein, 
one also needs an external constraint to determine the 
magnitude of the potential, as Eqs. (24) through (35) 
are invariant under the scaling K 1 T — ► aK 1 aT 1 and 
{b),(b 2 ) — > (b)/a, (b 2 )/a 2 (within the limits of variation 
of b). The simplest such constraint is obtained from the 
knowledge of the redshift of a giant arc (if one is present 
in the field, which fortunately is not an exceptional sit- 
uation). As for the source redshift distribution, (b) and 
(b 2 ) depend on both the redshift of the lens and on the 
Cosmology besides the shape of the source redshift distri- 
bution, so that it seems likely that the combination of the 
application of the reconstruction procedure just outlined 
to a number of lenses at different (known) redshifts with 
other lensing methods depending on the characteristics of 
the redshift distribution and constraining the Cosmology 
(see, e.g., Fort et al. 1997, and Bernardeau et al. 1997) 
will finally yield both the source redshift distribution and 
the cosmological parameters f2 and A. 

Although it seems less likely to occur in practice, let 
us now also consider the other possibility, namely that the 
critical zones are extended enough in the image plane so 
that they cannot be ignored in the reconstruction proce- 
dure. Comments made above and applying to this case 
(with appropriate transpositions) are not repeated here. 
It is likely that the inversion is still possible in princi- 
ple, although the reconstructed surface density becomes 
Cosmology dependent. For example, one can assume an a 
priori Cosmology and redshift distribution, which fixes the 
AB n factors in Eq. (32). Then Eq. (32) depends only on 
the lens mass distribution, which can be recovered in much 
the same way as described above along with (b) and {b 2 ) 
(if needed) which can be considered as free parameters for 
the mass distribution reconstruction. Then Eq. (32) can 
be solved for the redshift distribution for a given Cosmol- 
ogy and the previously found surface density and possibly 
the values of (b) and (b 2 ), and the procedure iterated until 
convergence is obtained. The weak points of this approach 
are that the limited size of the concerned regions in the 
image plane, the possible breakdown of the approximation 
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underlying Eq. (10) for some sources not individually iden- 
tified, the nonvanishing size of the background sources, 
the limited background sources number density and other 
observational sources of noise and bias are likely to make 
such a reconstruction procedure much more difficult to im- 
plement than the preceding one (not mentioning possible 
problems in the reconstruction of p{z$) from the integral 
equation). The strong point is that a single well-chosen 
lens, combined with the methods mentioned above to con- 
strain the Cosmology could in principle yield at the same 
time the cosmological parameters and the source redshift 
distribution, a rather tantalizing possibility. 

Acknowledgements. I thank Yannick Mellier and Peter Schnei- 
der for informal but informative discussions on the content of 
the paper. 
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